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We study as function of energy strangeness created in relativistic heavy ion collisions. We consider 
statistical hadronization with chemical freeze-out in both equilibrium and nonequilibrium. We 
obtain strangeness per baryon and per entropy in the energy range 8.75 ^ -y/ijvjv ^ 200^1 GeV. A 
baryon density independent evaluation of the kaon to pion ratio is presented. 

PACS numbers: 12.38.Mh,24.10.Pa,25.75.-q 



Production of strange hadrons in heavy ion collisions 
has been predicted to be sensitive to deconfinement. 
Firstly, there is the establishment by way of gluon fu- 
sion reaction gg — > ss of an abundant supply of strange 
quarks and antiquarks. Once the quark-gluon plasma 
(QGP) cools to the point of hadronization, in a second 
and independent step the final state hadrons are made 
in a recombination-fragmentation process Q. A spe- 
cific deconfinement feature is the enhanced production 
of strange antibaryons, increasing with strangeness con- 
tent, a feature seen in recent experiments 0- 

We analyze here experimental results in search for 
discontinuities in excitation of strangeness as function 
of reaction energy. The production yields of (strange) 
hadrons are studied in several experiments at the RHIC 
and at the SPS. Understanding of hadron yields in terms 
of phase space densities allows us to evaluate the global 
properties of all particles produced . The particle pro- 
duction can be well described in a very large range of 
yields solely by evaluating the accessible phase space size, 
when including many hadron resonances 0. 

In this work, we consider chemical equilibrium and 
non-equilibrium, i.e., we allow quark pair phase space oc- 
cupancies, for light quarks jq ^ 1, and/or strange quarks 
7s 7^ 1 iSl f-iid we require balance of strange and anti- 
strange quark content Q]. There are two independent 
fit parameters when we assume complete chemical equi- 
librium, the chemical freeze-out temperature T and 
the baryochemical potential (or equivalently, the quark 
fugacity Xq = e^^^^'^). Adding the possibility that the 
number of strange quark pairs is not in chemical equilib- 
rium, 7s ^ 1, we have 3 parameters, and allowing also 
that light quark pair number is not in chemical equilib- 
rium, we have 4 parameters. These three alternatives will 
be coded as open triangles (green online), open squares 
(violet) and filled squares (red), respectively, in the re- 
sults we present graphically. 

Statistical hadronization cannot be modeled com- 
pletely today, as we neither know all hadron resonances. 



nor do we know the required branching ratios of reso- 
nance decays. This introduces arbitrariness in the model 
which can lead to discrepancies between hadronization 
analysis results. To estimate the systematic error, we 
have performed a study varying the pion yield artificially 
by a factor 0.8 < f-w < 1.2. One can infer from this study 
that, if a relatively large hadronization temperature is 
reported, the presumption must be made that the sta- 
tistical hadronization program used does not produce as 
many pions as required by an extrapolation of resonance 
mass spectrum and resonance decay pattern. Our com- 
puted yields have been cross-checked (for the chemical 
equilibrium variant) as noted in acknowledgments. 

At RHIC, the baryon yield is found to contain many 
strange baryons. Not all results have been corrected for 
ensuing weak decays. We assume in our approach that 
50% of weak decays from S to A and from to S are 
inadvertently included in the yields, when these had not 
been corrected for such decays. We further assume that 
pions from such weak decays are not included in the ex- 
perimental yields, as these pions can clearly be shown to 
originate outside the interaction vertex. 

We have carried out a RHIC-200 (i.e., ^/smn = 
200 A GeV) Au-Au reaction analysis based on reported 
BRAHMS T,'8^, PHENIX ij, PHOBO^IO], STAR [ll 
iUda], results for tt, h'', p, p, K, K, K* , 0, 
yields. We have further reanalyzed the extensive RHIC- 
130 Au"Au results (compare both in order to ac- 
count for latest resonance yields and to make sure that 
there is no significant change introduced by the refine- 
ments made in the hadronization program. For the SPS, 
we take results obtained with Pb beams reacting with 
Pb stationary target at ^/snn = 8.75, 12.25, 17.2 GeV 
(projectile energy 40, 80, and 158 A GeV). We use here 
the SPS NA49-experiment 47r particle multiplicity re- 
sults Ullll, which include tt^, K, K, A, A, at 40, 
80, 158A GeV. We also fit (relative) yields of S, S, n, ft 
when available. Since we fit 47r-particle yields, no infor- 
mation about the collective fiow velocity is obtained. 
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TABLE I: The chemical freeze-out statistical parameters 
found for nonequihbrium (left) and semi equilibrium (right) 
fits to RHIC results. We show y^sjvjv, the temperature T, 
baryochemical potential ^t, strangeness chemical potential 
fj,s, the quark occupancy parameters 7, and 7s /iq, and the 
statistical significance of the fit. The star (*) indicates that 
there is an upper limit on the value of 7^ < e""/'^ (on left), 
and/or that the value is set (on right). 
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21.5 ±1 29.2 ±1.5 
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1.2 ±0.15 1.3 ±0.1 


1* 1* 
1.0 ±0.1 1.13 ±0.06 
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2.9/6 15.8/24 


4.5/7 32.2/25 



We present the RHIC freeze-out statistical parameters 
in table m These results are incorporating statistical and 
systematic errors in the data, along with errors in the 
statistical hadronization theory, arising from the above 
described uncertainty of the pion yield. The bottom line 
in table Q presents the statistical significance. The in- 
troduction of the full chemical nonequihbrium (on left) 
reduces by factor two the value of ■ Even when allow- 
ing for all systematic uncertainties and theoretical un- 
certainty regarding the pion yield, the nonequihbrium fit 
is a more compelling considering the (presently) more 
data rich RHIC-130 system. The confidence level of the 
chemical equilibrium approach remains here at 15%. At 
RHIC-200, the current data can be interpreted within a 
chemical equilibrium model, since the decisive results on 
relative yields of ( mult i) strange (anti)hyperons are not 
yet available. 

The statistical hadronization is expected to describe 
particle production well in presence of a sudden QGP 
breakup. Statistical hadronization approach is not ap- 
propriate when there is a lot of hadron-hadron rescat- 
tering in the final state, which may be the case at low 
SPS energies and below. In this case, kinetic models need 
to be applied, introducing a multitude of freeze-out con- 
ditions depending on the nature of particle considered. 
Accordingly, when we fit the low energy SPS results, a 
less favorable is found than at RHIC, or at the top SPS 
energy. However, the stability of the physical properties 
we extract when we subject the system to perturbations 
regarding pion yield indicate that results concerning the 
physical properties of the hadron source are reliable. 

With the statistical parameters fixed, the properties 
of the fireball can be computed, evaluating the contribu- 
tions made by each of the hadronic particle species. We 
first show how this works considering the energy stop- 
ping in figure ^ This result allows us to represent be- 
low our further findings as function of Ef^j^, the specific 
per baryon pair thermal energy available at the time of 
hadronization. This dependence substitutes for the de- 



0.6 



0.5 



0.4 



0.3 



A' 



A 



'U 



,A T, 

_l I I ! I I r 



Fitted parameters <^ □ t' x^' Z^' 

" ' q' ' s 



10 



100 



^ [GeV] 



FIG. 1: Fraction of energy stopping at SPS and RHIC; results 
are shown for 40, 80, 158A GeV Pb-Pb, 200A GeV S-W/Pb 
reactions and at RHIC for 65±65A GeV Au-Au interactions. 
Connecting lines guide the eye (color online). 



pendence on ^/sNN, the initial energy per baryon pair 
brought into the reaction. We see in figure that the 
chemical equilibrium fit (open triangles) shows counter- 
intuitive behavior. 

We note that the fraction of the energy per baryon pair 
which is initially thermalized is obtained from this result 
by adding the kinetic energy of the collective matter flow. 
Using at RHIC-200 as the average transverse flow veloc- 
ity {v'j_) — 0.45^ (and smaller values at smaller collision 
energies) we find a ~ 20% correction. Thus, we see that 
50 ± 5% of the energy per baryon carried into the reac- 
tion is initially made available to the thermal degrees of 
freedom, and this value is independent of the collision en- 
ergy in the entire SPS and RHIC range 9 < ^/snn < 200 
GeV. This means that the baryon stopping is two times 
larger then energy stopping. 

We are now ready to consider the final state specific per 
baryon strangeness yield. We evaluate strangeness yield 
s, the number of produced strange quark pairs and divide 
it by the similarly computed thermal fireball baryon num- 
ber content b. b is obtained summing the net (particle 
minus antiparticle) yields of all nucleons, (multi strange) 
hyperons and their resonances. The baryon number b is 
naturally conserved. Strangeness is predominantly pro- 
duced in the early stage of the reaction, when the density 
and temperature are highest, and there is little change 
in this yield during the late fireball evolution. Thus, s/b 
ratio probes directly the extreme initial conditions. Con- 
sideration of this yield ratio eliminates the absolute yield 
normalization parameter (sometimes but not always 're- 
action' volume), as well as some uncertainties originating 
in the experimental particle yield data, which propagate 
through the analysis of experimental data. 

In figure 12 we show s /b as function of the thermal 
specific energy content Ef^j^. Strangeness yield is con- 
tinuous in the entire energy domain as is shown by the 
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FIG. 2; Strangeness per thermal baryon participant, s/b as FIG. 3: Strangeness per entropy, s/S, as function of the 



function of thermal specific energy content j 



thermal specific energy content E^j 



solid line drawn to guide the eye. The rise of specific 
strangeness yield is slightly faster than linear with energy, 
as is indicated in the figure insert. (The appearance of an 
exponential shape is due to logarithmic energy scale in 
figure|2l. The results for SPS energy range are in quanti- 
tative agreement with predictions made assuming a QGP 
state of matter and gluon fusion strangeness formation 
mechanism, compare figure 38 in Ref. [23|- The SPS 
specific strangeness yield extrapolates smoothly to the 
RHIC energy range. The reaction mechanism producing 
strangeness and stopping baryon number are evolving in 
parallel yielding a smooth change in the ratio of both 
variables. 

The highest energy ^/snn=200 GeV data point in fig- 
ure El seems to be slightly lower than the presented ex- 
trapolation predicts. We think that this yield will in- 
crease once we have included in the statistical hadroniza- 
tion analysis the hyperon yield ratios S/A,A/p. The 
presence, in the global fit, of these particle ratios will 
increase the fitted strangeness yield, provided that their 
measured values are comparable in magnitude to those 
reported at ysjv]v=130 GeV. 

An important feature of the RHIC experimental re- 
sults, shown in figure HI is the large strangeness yield 
per participating baryon. Although it has been early on 
noted that at RHIC- 130, up to 8 strange quark pairs per 
baryon are produced [13, little attention has been given 
to this high yield. At RHIC-2000 each interacting baryon 
pair produces about 20 strange quark-antiquark pairs. 

An interesting question is how this great increase in 
yield compares to the production of particles in general. 
The Wroblewski ratio Ws = 2{ss) / (uu + dd) is often used 
in such a comparison. Recently it has been realized that 
this ratio can be artificially enhanced by high baryon 
density which can suppress light quark pair yields |l8| . 
Therefore we here compare the strangeness production 
to the global entropy S yield. We evaluate S in the 
same way as we have obtained other global properties. 
There is an additional nonequilibrium entropy term to 



be allowed for when chemical non-equilibrium prevails. 
Both entropy S and strangeness s are nearly conserved 
in the hydrodynamical expansion of QGP and increase 
only moderately in the hadronization of QGP. The ob- 
served ratio s/S is established by microscopic reactions 
operational in the early stages of the heavy ion collision. 

Strangeness per entropy, s/S, is presented in figure O 
as function of the specific thermal energy content. The 
horizontal line is the maximum SPS yield base. We note 
the modest smooth rise in the SPS energy domain. Most 
interestingly, at RHIC, as compared to SPS, an unex- 
pected 50% increase in s/S is noted, allowing for chemi- 
cal nonequilibrium. The excitation of strangeness seems 
to rise faster with energy than the production of entropy. 
A possible explanation of this phenomenon is that the 
hot initial state, in which the threshold in energy for 
strangeness formation has been overcome, lives longer 
when formed at RHIC conditions. However, it appears 
important to fill the energy range between SPS and RHIC 
with data in search of a new physics energy threshold. 

Based on the study of the K'^/tt^ particle yield ra- 
tio 1^, such a new physics energy threshold had been 
expected below 40 A GeV, at the lowest SPS energy, i.e., 
below the energy range we could consider in the global 
hadron freeze-out analysis. Inspecting available particle 
yields, we note that their dependence on the unknown 
baryon density is significant, mimicking new physics. 
This baryon density effect can be greatly reduced con- 
sidering the product of particle and antiparticle yields. 
In this case, the chemical potential cofactors in particle 
yields cancel, being inverse of each other, and the baryon 
density effect largely disappears. To study kaon to pion 
ratio we thus form: K/n = yj [K+K^)/{t:+tt~), and 
effectively 'remove' baryon density effect present in the 
individual ratios /tt^ and /tt^ . Using the NA49 
data set jl^, we present the K/tt ratio in figure 0] The 
filled squares are for the 4-7r full phase space data set. 
The filled circles are for the central rapidity results, in- 
cluding RHIC-130 and RHIC-200 data The central 
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FIG. 4: K/tt ratio as function of collision energy. Filled 
symbols are for nuclear and open for K'^/n'^ measured in 
elementary pp collisions. Squares denote the full multiplicity 
ratio and circles the central rapidity yield ratio. 



rapidity NA49 results are from figure 6 in The pp 
charged A'+/7r+ > K/tt background is shown by open 
squares, and we indicate for = 1800 the p-p TEVA- 
TRON point which is at an energy beyond the range 
considered. 

The lines, in figure 0] guide the eye to the trends of 
the alternate gradient synchrotron (AGS) and SPS re- 
sults. These trends of behavior intersect within the SPS 
low energy range below the lowest NA49 SPS point which 
was obtained at SOAGeV. However, this appears here to 
be a smooth transition from a rise to saturation of the 
K/tt production. There is a clear enhancement of the 
K/tt ratio at RHIC compared both with the AGS/SPS 



trend and with the /tt^ > K/tt measured in elemen- 
tary pp collisions. This enhancement coincides with the 
specific strangeness per entropy, s/S enhancement seen 
in figure Ul 

We have shown that the strangeness per baryon exci- 
tation in the fireball of dense matter formed at SPS and 
RHIC energies is a smooth function of energy, but we find 
a step-up in strangeness per entropy between the SPS and 
RHIC energy ranges, also visible in the rise of the K/tt ra- 
tio. Thus there are two energy domains to investigate for 
threshold behavior, the 35-1-35 GeV RHIC range and the 
20A GeV on fixed target low SPS energy, where the K/tt 
ratio saturates. We further noted the quantitative agree- 
ment in the strangeness yield with predictions made for 
the SPS energy range, assuming QGP production mecha- 
nisms. We have shown that 20 strange-antistrange quark 
pairs are made per colliding baryon pair at the top RHIC 
energy. 

Note added More technical details are now available in 
Ref. |22|- A systematic study of SPS results by a yet 
another group has just appeared and the numerical 
results shown there agree with our SPS results. Fig 10 
shows that 7^ = 1 is actually a local fit maximum for the 
top SPS energy. 
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